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1 Introduction 



The decays r]^'^ tt+tt^/^/^ are interesting in several respects. First, they involve contributions 
from the box-anomaly of quantum chromodynamics. Second, they probe the transition form 
factors of the 1] and 1]'. In principle, the decays are suited to test whether double vector meson 
dominance is indeed realized in nature, which is also an important issue for the anomalous 
magnetic moment of the muon and kaon decays [1]. Moreover, since the r]' is closely related 
to the axial U(l) anomaly of the strong interactions, one can study the phenomenological 
implications of the anomaly at low energies. 

On the experimental side, there is renewed interest in t], t]' decays which are investigated at 
WASAOCOSY [2], MAMI [3], KLOE [4,5] and by the VES collaboration [6,7]. There is thus 
the necessity to provide a consistent and uniform theoretical description for these decays. 

In this respect, the combination of the chiral effective Lagrangian which incorporates the 
symmetries and symmetry-breaking patterns of QCD in combination with a coupled-channels 
Bethe-Salpeter equation (BSE) that takes into account final-state interactions in the decays 
and satisfies exact two-body unitarity has been proven very useful. In a series of papers, this 
approach has been successfully applied to the hadronic decay modes of t] and rj' [8-10], and the 
anomalous decays r]^'^ — > [11] and r]^'^ tc^tt^'j [12]. 

Of particular interest is the last work [12] which we extend here to off-shell photons since the 
process i]^'"* TT^n~l^l~ can be regarded as the two-step process r]^'^ n^n^^* n^n^l^l^. 
It is worthwhile mentioning that the conventional vector dominance picture with energy- 
dependent widths in the vector meson propagators can be shown to be in contradiction to 
the one-loop result of chiral perturbation theory (ChPT) [13], the effective field theory of the 
strong interactions. The present approach, on the other hand, satisfies theoretical constraints 
such as anomalous Ward identities, electromagnetic gauge invariance, exact two-body unitarity 
and matches in the low-energy limit to one-loop ChPT. Resonances are not taken into ac- 
count explicitly, but are rather generated dynamically through the iteration of meson-meson 
interactions. 

This work is organized as follows. In the next section we present the general structure of the 
amplitude, while in Sec. O the one-loop result of these decays within ChPT is derived. Some 
details of the chiral unitary approach are presented in Sec. H] and the results are discussed in 
Sec. [5l We summarize our findings in Sec. El The full list of relevant 0{p^) counter terms is 
relegated to the appendix. 

2 General structure of the amplitude 

The decays rj^'^ tt^tt^I^I^ (l^ represents either or /i^) are depicted in Fig. [H where we 
also introduce the four-momenta of the particles. The invariant matrix element of the decay 
has the generic form 

= -tee^'"'^k,ptp-A{s+^,s^„s.,)^u{q~,a){-zeYHq^,cr') , (1) 

with spin indices a, a' and 74(s+_, summarizing all contributions to i]^'^ — > n^rr^'y* 

(represented by the blob in Fig. [1]). The Mandelstam variables are defined as 
follows: 

s+_ = (p+ + , = (p+ + kf , s.^ = {p- + kf . (2) 
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Figure 1: General structure of the process r]^'\P) — > 7r^(p'^)7r~(p~)/^(g^)/~(g~). The blob 
symbolizes the amplitude A{ri^'^ 7r+7r^7*). The four-momentum of the intermediate photon 
is denoted hy k — P — — = + q~ with k'^ > 0. 



As a consequence of C invariance A{s+^, is symmetric under the exchange <-> 

Clearly, the decay i]^'^ 7r^7r~l~^l~ proceeds via the two-step mechanism r]^'^ — >■ 71~^tt~^* 
followed by 7* . Defining (in accordance with [14]) the n-body phase space element 

dK{P;Pi, ...,Pn)= s^'^ (p - J2p?) n (^^5% (3) 

and making use of the factorization 

d$4(P; g+, q-,P^,p-) = d<^s{P; k,p+,p-) d^2{k] g+, q-) {27rf dk^ (4) 

one finds the following relation between the differential decay width of 77*^'^ — > tt+tt^Z+Z" and 
the differential widths of the two sub-processes 77*^'^ 7r+7r'~7* and 7* see e.g. [15]: 

dT{r]^'^ ^ TT+TT-l+r) = dr{ri^'^ ^ 7r+7r-7*) dT{Y ^ - ^= dk^ . (5) 

After integration over the dilepton phase space (PS//) one arrives at 

/ dr{r]^'^ n+n-l+r) = dr{r]^'^ 77+77-7*) r(7* ^ /+/-) - dk'' (6) 

with 



X Oi r—rf 2mf\ / 4m? , , 

r(,.^,*r) = -v^(i + ^)^i-^ . „ = (7) 

The task of the current work is to calculate yl(s_|__, s_-y) (or, equivalently, the amplitude 
A{rj^'^ — > 7r"'"7r-7*)) within a chiral unitary approach. 



3 One-loop calculation 

In this section we present the result of the full one-loop calculation of the amplitude for 77*^') 
7r+7r~7* in U(3) ChPT generalizing the one-loop result of [12] for the decay amplitude 77*^'^ 
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7r+7r~7- Here we will restrict ourselves to compiling the necessary formulae and outlining the 
basic steps of the calculation. For details we refer the reader to [12]. 

The amplitude A{r]^'^ 7r+7r~7*) involves the totally antisymmetric tensor e^"^"^ and is 
thus of unnatural parity. At leading chiral order, the pure SU(3) process 1]^ 7r+7r^7* is 
determined by the chiral anomaly of the underlying QCD Lagrangian. Within ChPT the chiral 
QCD anomalies are accounted for by the Wess-Zumino-Witten (WZW) action [16-20] 

S^zw = j d'x e'^^'^iU^d^UU^djUU^dkUU^diUU^dmU) (8) 



Ms 

e 
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where we have displayed only the pieces of the action relevant for the present calculation. The 
octet of Goldstone bosons (vr, K, tjq) and the singlet field rjo are collected in the matrix valued 
field which enters into U = exp{zv^0//}, where / is the pseudoscalar decay constant in the 
chiral limit. The expression (. . .) denotes the trace in flavor space, is the photon field, and 
Q = |diag(2, —1, —1) represents the charge matrix of the light quarks. The integration in the 
first line of Eq. ([H]) spans over a five- dimensional manifold M5, whose boundary is Minkowskian 
space, and the U fields in this integral are functions on M5. The additional fifth coordinate is 
defined to be timelike and the convention for the totally antisymmetric tensor is e0^234 = _|_x^ 
see [17, 19, 20] for further details. 

The inclusion of the singlet field rjQ and, consequently, the extension of SU(3) ChPT to the 
U(3) framework introduces additional, non-anomalous terms of unnatural parity at chiral order 
0{p^). The only term relevant for this work at O(p^) reads 

4? = -te e^^^^d^A, {dJJdpU^Q + d^U^d^UQ) , (9) 

where VF3 is a function of r^o, ^3(^70//); which can be expanded in the singlet field with 
coefficients w^^ that are not fixed by chiral symmetry. Parity conservation implies that is 
an odd function of r/o- 

In addition to the leading-order tree level contributions derived from Eqs. ([H]) and Q there 
are next-to-leading order chiral corrections from one-loop graphs, decay constants, rj-rj' mixing, 
and wave function renormalization which involve terms both from the 0{p^) + 0{p'^) Lagrangian 
and the C(p^) Lagrangian of natural parity with couplings v^^ '' and j3'f\ respectively. The full 
list of terms up to 0{p^) can be found, e.g., in [21]. Finally, the process rj^'^ 7r+7r~7* receives 
contributions from counter terms of the unnatural parity 0{p^) Lagrangian, which also absorb 
the divergences of the one-loop integrals. 

Fig. [2] shows the pertinent one-loop diagrams contributing to ?7*-'^ — 7r"'"7r~7* (except for 
contributions from wave function renormalization). The full one- loop result reads 

A('-'-^\v^'^ - 7r+vr-7*) = -eA:,6.p>,-6^-^ ^^^^^^^^^^^ ^^^-^"^^ , (10) 
where is the polarization vector of the virtual photon and the coefficients (3^J^~l°°^^ are given 
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Figure 2: One loop diagrams contributing to the process 77^'^ — > 7r+7r~7*. The empty squares 
denote vertices from the 0{p^) Lagrangian of unnatural parity, whereas vertices from the leading 
order Lagrangian of natural parity are indicated by a filled circle. 
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- 3A(m^) - -A(m^) + /(m^, m^, s+_) + -/(m^, 



(11) 



In this expression we have perturbatively substituted the pseudoscalar decay constant in the 
chiral Umit, /, by the physical decay constants F^^, F^, F^r of tt, 77, rj', respectively, and employed 
the abbreviations 



_ 1 f2 1 



'3 > 



P5,18 — P5 + 2/^: 



18 



(12) 



The loop integrals are calculated using dimensional regularization and the pertinent regular- 
ization scale is denoted by /i. The finite parts of the loop integrals are given by 



A{m^) 



dH 



2 2 
m m 



(13) 



and 



I{m^,m'^,p^) =g^{ - (P^ - (m - mf){p'^ - {m + mf)Gmm{p^) 

+ (p'^ + m'^ - m^)A(m^) + (p^ - + m^)A(m^)} 



(14) 



+ 



1447r2 



(p — 3m — 3m ) , 
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where Gmfn is the finite part of the scalar one-loop integral 

■2-. . . d'^l 



(2txY {P -m? + ie){{l - - m? + ie) 
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— artanh ■ 



finite 



[m + m)2 — 



(15) 



Amm(p^) = {{m - mf - p^) ((m + mf - p^) . 
The integral /' is defined via the subtraction 

/'(m^ m^, p^) = I{m^, m^,p^) - /(O, 0) 



(16) 



which guarantees chiral power counting for loops involving the rj' . Since the mass of this heavy 
degree of freedom does not vanish in the chiral limit, its presence can in principle spoil the 
chiral counting scheme. However, it has been shown in [12] that all power-counting violating 
contributions to the process rj^'^ 7r+7r~7* can be absorbed into a redefinition of the low-energy 
constant w^''] the renormalized value is denoted by w^''^ . 



The last terms in the expressions for 13 



(l-loop) 

,(') 



in Eq. ( fTTj) summarize the contributions of 



counter terms from the 0{p^) Lagrangian of unnatural parity. The relations between the 

tZ;^(,^) and the numerous couplings of the Lagrangian of sixth chiral order 
are given in the appendix. 



constants w^^}^, w^"^^ 



(')' 



4 Chiral unitary approach 

From the analysis of various rj and 1]' decays, see e.g. [9, 11, 12], it has become clear that 
resonances and unitarity corrections due to final-state interactions are a necessary ingredient 
for the realistic description of these processes. One example is the pronounced peak structure 
caused by the p(770) resonance in the tt+tt^ spectrum of rj' 7r+7r^7 [22,23] (see also Fig. E]). 
Hence, a conventional loop-wise expansion within ChPT is usually not sufficient to successfully 
describe rj and — in particular — rj' decays. 

Instead of taking resonances into account explicitly, as e.g. in [15,24,25], we prefer to work 
within a chiral unitary approach which combines ChPT and a non-perturbative resummation 
based on the Bethe-Salpeter equation (BSE). In this framework the resulting multi-channel T- 
matrix of meson-meson scattering satisfies exact two-body unitarity. Our approach has the fur- 
ther advantages that electromagnetic gauge invariance is automatically maintained, anomalous 
chiral Ward identities are satisfied, and the result matches to one-loop ChPT in the low-energy 
limit. Resonances are generated dynamically and are identified with poles of the T-matrix in 
the complex energy plane. 

Since this approach has already been discussed in detail in [11,12] we will only recapitulate 
the basic formulae here. From the effective Lagrangian up to fourth chiral order one extracts 
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the partial wave interaction kernel Ai for meson-meson scattering which is then iterated in the 
BSE 

Te = Ae-AeGTe. (17) 
The diagonal matrix G collects the modified scalar loop integrals 



(18) 



where we have added a subtraction constant ammifJ') to the integral Gmm defined in Eq. ffT^ 
which varies with the scale /i in such a way that Gmm is scale-independent [26]. After adjusting 
the occurring parameters the partial-wave T-matrix resulting from the BSE accurately describes 
the experimental phase shifts in both the s- and p-wave channels [11,27]. 

The implementation of non-perturbative meson-meson rescattering generated by the BSE 
in the amplitude A{r]^'^ 7r+7r^7*) is accomplished in the same way as in [12]. The pertinent 
graphs are shown in Fig. [H] and the corresponding amplitude is added to the one-loop result 
presented in in the previous section. We point out that a possible double counting of one-loop 
contributions, which in principle arises since the diagrams (a) and (c) in Fig. [3] incorporate also 
one-loop terms, has been properly taken care of. The amplitude corresponding to the diagram 
in Fig. [3^ is given by 
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Y^l^^T^''' I{ml ml C^) fi,^-^\s^.) (19) 



with 



1: 



(CCa),K°K° 
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4^2, 



K 



.,(2) 
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V' 



1. 
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1 
6 



(CCa),K^ 



VQ — 4:bTT~W 



2„.,{l)r 
3 



-2W^(m 



+ 



4 



4^6 
3 ^ 

/^5,18 



„,(2) 



18 



F2 



„,(2) 



(20) 



The symbol Yl' Eq. ffTI?|) denotes summation over the meson pairs n'^n , K^K and K^K^ 
and Tp°'~'^^ represents the p-wave part of T-matrix for scattering of a meson pair a into a meson 
pair b. It differs from the solution of the BSE, Tp, only by kinematical factors, see [11] for 
details. The loop integral / is given by 



I{m'^,m^,p'^,Gmfh) = I {m'^ , , p^) + (p^ — 3m^ — 3m^)C„ 
with / defined in Eq. ( fT^ . In order to keep the notation compact we set 



Ctt = G, 



Gk = G, 



n — n 



G-jxrji — Cn 



(21) 



(22) 
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(b) (d) 7%k 



Figure 3: Set of meson-meson rescattering diagrams which contribute to the process r]^'^ — >■ 
TT+Tf-^* and are taken into account in this approach. The empty squares denote vertices from 
the 0{p'^) Lagrangian of unnatural parity whereas vertices from the leading order Lagrangian 
of natural parity are indicated by a filled circle. 
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Note that the definition of / shghtly differs from that of Ji in [12], where the constant C was 
chosen to be the coefficient of instead of {p^ — "im? — 3m^). Here, we prefer to work with the 
decomposition in Eq. fl?Il) since then the regularization scale dependence of / can be completely 
absorbed into the constant C. We point out that in an effective field theory framework one is 
free to arbitrarily modify the analytic piece of an amplitude by adjusting unconstrained counter 
terms. 

The diagram in Fig. [3b produces the amplitude 

A(cc^^^i') = -eKe.vlp-,^^'''' ^YH""'^" i{ml,ml,k\C^) 

X f^'^--^\e) i{mlmlk\C^) +fir-''"\e) i{m^^^^^ , (23) 

where the coefficients 7 ''('^f are given by 



^(CCb),^± ^ (CCb),n± ^ „ ^(CCb),K± ^ _^(CCb),KOKO ^ 

iTj ir/' ^ 1 iTj Iri 2 ' 



(CCb),K± (CCb),K^K° t 2 2\ { a(^5,18 i 
%' =-%' = ^i'TUj, - m^) Uy^^ (24) 



The amplitude corresponding to graph (c) in Fig. [3] is given by 

1 



+ i\ml,ml„s,„C^,,) ff''''-'^'-'\s_,)]y (25) 

where the integral /' is defined analogously to /', Eq. ( |T6l) . by 

/'(m^m^p^C*™^) = /(m^m^p^C^^) - J(0,m^O,C'm^) . (26) 

Finally, we include the diagram with two insertions of iterated meson-meson rescattering. 
Fig. [3li. The corresponding amplitude reads 



[f^r-^\e) i{mlmlk\C^)+f^'-^^\k') i{ml,ml,e,CK)] (27) 



X 
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with coefficients -y'^^'"'^^'"'* symmetric under the interchange a -^h 



(2CC),'K^,K± _ _^(2CC),TT^,K°K° _ _\ (2CC),K^,K°K° . 



/?7 /r; n h 



(2CC),-K^,K^ _ _ (2CC),n'^,K°K° _ _1 (2CC),K'^,K°K° 
Iri 7?j' 2 7t;' 

and zero otherwise. 

5 Results 

The chiral unitary approach discussed in this work involves several free parameters which must 
be fixed from experiment. On the one hand, there are the coupling constants of the chiral 
Lagrangian which can be grouped into coefficients of the natural parity part of OijP^ + Oij?^ 
and Oijy")^ v'p and (i^p ^ respectively, and coefficients of the unnatural parity part of C(p^) and 
C(p^), wf"^ and wf \ respectively. On the other hand, there are the subtraction constants a and 
C in the loop integrals whose values correspond to a specific choice of the infinitely many higher 
order counter terms neglected in this non-perturbative approach. For consistency with previous 
work [9] the coupling constants of the Lagrangian of natural parity and the subtraction constant 
a^^T^'^~^^ in the isospin one p-wave tttt channel are fixed by a fit to the hadronic decay modes 
of T] and 77', r\^'^ — > Stt and Vj — > t^tttt, and the phase shifts of meson-meson scattering. This fit 
is in very good agreement with the bulk of the available experimental data. The subtraction 
constants in the other meson-meson channels do not have any relevant impact on the discussed 
data and can be set to zero for our purposes. The pseudoscalar decay constants are set to 
= 92.4 McV, = 1.3F^, and F^/ = I.IF^ [11,12]. 

The couplings of the unnatural parity part of the Lagrangian and the subtraction constants 
C are taken as free parameters, which are constrained by fitting to the available spectra and 
widths of the decays r/*^'^ — > 7r"'"7r~7. It turns out, however, that in order to achieve agreement 
with the experimental data, only the subtraction constant in the pion loops, Cjr, is required to 
have a non-vanishing value, and we set all other subtraction constants to zero for simplicity. 
To further reduce the number of parameters and for consistency with previous investigations 
[11,12], we also set the renormalized coupling constant w^^^ of the unnatural parity Lagrangian 
at O(p^) to zero. We have confirmed that small variations in ^3^^^ do not alter our conclusions. 
Finally, the combinations of 0{p^) unnatural parity couplings denoted by which do not 

contribute to processes with on-shell photons and thus cannot be constrained by 77^'^ — > 7r"'"7r~7, 
will be neglected for the time being. Changes of the results due to non-zero values of these 
coefficients will be discussed at the end of this section. To summarize, there are five parameters, 
Ctt, wlj^\ wij'\ 'w^\ and w^f'' which are constrained by fitting the decays rj^'^ — > 7r'*"7r~7. 
Afterwards, we can predict the spectra and widths of 77*^'^ — > TT'^n~l'^l~ within this approach. 

The data of 77^'^ — > n'^n'^ involve the partial decay widths [14] and the di-pion spectra 
from [22,23,28,29]. In order to perform a global least-squares fit to these different data sets we 
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employ the following definition for the x^-function: 



d.o.f. N{Y,i rii-p) ^ Hi 

where is the number of observables and p the number of free parameters in the approach. 
The quantity xl is the standard x^-value computed for the i-th data set with rii data points. 
The above definition was introduced in [30] to equally weight each data set and to prevent, 
e.g., sets with only one data point (such as decay widths) from being dominated by sets with 
many data points (such as spectra). 

In order to quantify an error for our analysis we employ the condition [14] 

X ^ Xmin I (30) 



d.o.f. d.o.f. d.o.f. 

where Ax^ is derived from the p-value of the probability distribution function. One finds 
that in the present investigation employing Ax^/d.o.f. = 1.08 corresponds to the la confidence 
region. Strictly speaking, this standard definition of a confidence region, Eq. ( 130|) . holds only 
if the fit is performed to just one observable and the fit function is linear in the fit parameters. 
Although both constraints are not fulfilled here, one can expect Eq. f l30l) to be a reasonable 
approximation in the vicinity of the minimum of the x^-function, see also [31]. We have signif- 
icantly improved our fitting routine compared to the previous investigation [12] and performed 
a large number of fits so that the la confidence region is populated by about 1000 qualitatively 
different fits providing a realistic estimate of the theoretical uncertainty within this approach. 

In Figs, m and the result of the calculation is compared to the available experimental 
spectra which are given in terms of the photon energy for rj —>■ 7r+7r~7 and in terms of the 
invariant mass of the tt+tt^ system for rj' 7r+7r^7. The solid line corresponds to the best 
fit with an overall x^/d-o.f. = 2.23, the error bands indicate the la confidence level. For 
the r]' decay the agreement with the two experimental spectra from [22, 23] is very good as 
already observed in [12]. The experimental situation for the t] decay is not as consistent as 
for 1]' 7r+7r^7. First, both the spectra published in [28,29] have not been corrected for the 
detection efficiency which is given separately in [28], but must be deduced in [29]. Also, in 
both experiments it is impossible to quantify the systematic error resulting from the correction 
of the detection efficiency which introduces an uncontrolled uncertainty in the data. Second, 
when taking into account the two data sets from [28] and [29] simultaneously in the fit, it turns 
out that they are not fully consistent, at least without knowledge of the complete systematic 
errors. As a consequence, the major part of the total x^/d.o.f. value is due to the disagreement 
between the two data sets. In fact, the best fit (solid line in Fig. H]) must be considered as 
a compromise of [28] and [29], so that under these circumstances a total x^/d.o.f. close to 1 
cannot be achieved. If, however, only one of the two spectra is included in the fit, a total 
X^/d.o.f. ~ 1 can be obtained. In this context, further experimental investigations — such 
as [2] — with substantially improved accuracy should lead to a more consistent picture of the 
rj —>■ 7r+7r^7 spectrum. 

The numerical results for the branching ratios and decay widths of r^'-'-' 7r+7r^7, 7]^''^ — *■ 
7r~^7r~l~^l~ are shown in Table [TJ The central values of our results correspond to the fit with 
minimal x^? the error bars reflect the la confidence region given within our approach. The 
agreement with the decay modes involving on-shell photons, which have been taken as input to 
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Figure 4: Photon spectrum of ?7 — 7r+7r^7 compared to experimental data from [28] (left) 
and [29] (right). The solid line corresponds to the fit with minimal x^, the error band indicates 
the la confidence region. For comparison with the experimental data points, the curves have 
been multiplied by the experimental detection efficiencies, hence the different shapes in the two 
plots. 
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Figure 5: Invariant mass spectrum of the 7r"^7r~ system in rj' —>■ 7r+7r^7 compared to experi- 
mental data from [22] (left) and [23] (right). The solid line corresponds to the fit with minimal 
X^, the error band indicates the la confidence region. All curves are normalized to the integral 
of the experimental histogram. 
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Table 1: Results for the branching ratios (BR) and widths (F) of the decay modes under 
consideration compared to experimental values and the theoretical analyses [24,25] and [15]. 
See text for further details. 
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-20.ll?^5^ 
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1-^-11.8 


0+18.5 

-^•^-23.8 


1 q+0.7 



Table 2: Numerical values of the fitted parameters at the regularization scale /x = 1 GeV. The 
central values correspond to the fit with minimal the error ranges are given by the la 
confidence region. 

the fit, is very good. The numerical values of the fit parameters, i.e. the counter terms w^J,, 

w^^l) and the subtraction constant C^, are compiled in Table [21 Having fixed all parameters 
from data, we can make predictions for the decays into vr+vr" and a lepton-antilepton pair. 
Up to now, the only branching ratio of this type which has been determined experimentally 
is ?7 ^ 7r+7r"e+e~. We compare our result with the very recent experiment [32] which has 
improved precision compared to the PDG number [14] and we observe nice agreement. 

Moreover, we can compare our results with those of [24,25] and [15]. In [24,25] a chiral 
Lagrangian with explicit vector mesons is used to calculate both the decay widths and spectra 
of ?7 — s> 71+71- --^ and 77 71+71-e+e- . As shown in Table [1] and Figs. El [7] the agreement 
with our results is very good. However, it should be remarked that the results presented 
in [25] depend sensitively on the numerical values employed for the meson masses. Using 
the final expression Eq. (4) in [25] and inserting up-to-date meson mass values from [14], 
one computes r{ri 71+71-e+e-) = 403 meV instead of 380meV as given in [25]. Also, the 
invariant mass spectra shown in Figs. El [7] are rescaled accordingly. In [15], on the other 
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Table 3: Branching ratios of the decay modes into tt+tt^/+/^ with respect to the tt+tt^7 decays. 
The experimental value quoted in the third column is taken from [14]. The relative accuracies 
of the theoretical results are given in the last column. 

hand, a meson exchange model has been employed to calculate numerous decay modes of light 
unfiavored mesons. Despite dissimilarities between [15] and our approach the numerical results 
are in reasonable agreement. We point out that — in contrast to our work — no theoretical error 
estimates are given in [15,24,25]. 

The ratios between the tt+tt^/+/^ and tt+tt^7 decay channels are given in Tab. [HI The small 
theoretical uncertainties for the decays into an e+e^ pair are further reduced down to about 
1% in these ratios, while the theoretical accuracies for the ^'^ ^~ decay ratios remain roughly 
unaffected. This indicates that the Tr"'"Tr~e"'"e~ and Tr"'"Tr~7 decays are correlated which can be 
traced back to the shape of the QED part r(7* — > e^e~) / {k'^\/l^) in Eq. ([6]) describing the 
transition 7* e+e~. This function possesses a pronounced peak at the virtual photon mass 
/Cg = (1 + \/21)ml ^ 5.6 ml and projects out the values of the subprocesses i]^'^ tt+tt~7* 
at /Cg — close to the photon on-shell point k"^ = 0. For the fj^~^fi~ decays, on the other hand, 
the respective value k'^ = {1 + \/21)m'j^ ^ 5.6 is relatively far apart from fc^ = so that 
these decays are not immediately correlated to the 7r"'"7r~7 decays. We observe that for photon 
virtualities which are not too close to the upper boundary of phase space the rate r(?7*^') — >• 
7r+TT~7*) in our approach can be very well approximated by a Gaussian of the form T{r]^'^ — >■ 
7T+TT-7) exp(-A;VA2) with A = (97.8l^|) MeV and A = (l67.3+^|) MeV for the t] and V decay, 
respectively. In combination with the sharply peaked QED part the dependence on the small 
variations in A is further reduced in the branching ratios T(r]^'^ — > Tr+Tr~e"'"e~)/r(77(') — > Tr"'"Tr~7) 
resulting in the small relative uncertainties of about 1% mentioned above. 

In Figs. and [7] we present our predictions for the tt+tt^ and invariant mass spec- 
tra, respectively. The lepton-antilepton spectra are strongly peaked right above threshold, so 
for illustrational purposes we have multiplied these spectra by a factor k"^ which reduces the 
otherwise extremely pronounced peak. Due to the tiny branching fractions of the decays into 
7r+TT^yU+/i~ it will be experimentally very challenging to measure these kinds of spectra. The 
spectra of the decays involving an electron-positron pair, however, are likely to be probed at 
the ongoing experiment [2] at COSY-Jiilich. 

We reconfirm the findings of [12] regarding the importance of the different coupled channels 
diagrams in Fig. [3l The by far largest contribution to the decay amplitude stems from tt+tt" 
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(MeV) ^ (MeV) 

Figure 6: Predicted invariant mass spectra of the tt+tt^ system in the different decay modes. 
The sohd hnes represent the fit with minimal the error bands indicate the la confidence 
region. The result of [25] is represented by the dashed line in the upper left plot. 
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Figure 7: Predicted invariant mass spectra of the lepton-antilepton pair in the different decay 
modes. The sohd hnes represent the fit with minimal x^? fhe error bands indicate the la 
confidence region. For illustrational reasons the spectra are multiphed by a factor k'^. The 
result of [25] is represented by the dashed line in the upper left plot. 
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Table 4: This table illustrates how the uncertainties of the results grow if variations of k'^- 
dependent counter terms are taken into account. 

final-state interactions, cf. Fig. [3K, whereas the diagram in Fig. [3li, which mimics the simulta- 
neous exchange of two vector mesons within our approach, yields only small corrections. This 
is in contrast to the assumption of complete vector meson dominance. 

The decay modes r]^'^ 71+71^1+1- involve 0{p^) counter terms which generate contributions 
proportional to k"^ and do not contribute to the decays with on-shell photons. Consequently, 
they cannot be fixed by fitting to r]^''^ 7r+TT~'j data. In order to examine their impact on the 
results for 7]^'^ — >■ Tr+Tr"/"*"/" we have varied their values in the range (—10 . . . + 10) x 10^ GeV~^ 
which is motivated by the size of the other 0{p^) couplings, cf. Table [21 The enlargement of 
the error ranges for the branching ratios and widths following from this variation is tabulated 
in Table m It turns out that the influence of the k"^ terms is rather mild for the decays involving 
an electron-positron pair owing to the fact that the spectra of such decay modes are strongly 
enhanced at small fc^, cf. Fig. [71 For the decays into 7r+n~n+^~ , on the other hand, where k"^ is 
bounded below by 4m^, the uncertainties from the la confidence regions are roughly doubled 

by taking into account the counter terms tv'"^}^ . 

Finally, we show in Tab.[nithe ratios r(?7*^'^ — ^ n+n-l+l-) /T{r]^'^ — n+ir-'j) in the presence 
of the k'^ terms. The relative uncertainties for the fj^+fi~ decays are again approximately doubled 
with respect to Tab. [3l if these counter terms are taken into account. 

6 Conclusions 

In this work we have investigated the decays rj, rj' — >■ tt+tt-I+I- within a chiral unitary approach 
based on the chiral effective Lagrangian and a coupled-channels Bethe-Salpeter equation. Uti- 
lization of the chiral effective Lagrangian guarantees that symmetries and symmetry-breaking 
patterns of the underlying theory QCD are incorporated in a model-independent fashion. In 
particular, contributions due to chiral anomalies enter through the Wess-Zumino-Witten La- 
grangian. Besides, counter terms of unnatural parity at leading and next-to-leading order are 
also taken into account. 
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Table 5: Branching ratios of decay modes into 7r^7r~l^l~ with respect to 7r"'"7r-7 decays when 
A;^-dependent counter terms are taken into account. The experimental value quoted in the third 
column is taken from [14]. The relative accuracies of the theoretical results are given in the 
last column. 

We have first performed a full one-loop calculation in ChPT. However, unitarity effects due 
to final-state interactions are important in rj and, in particular, in rj' decays and must be treated 
non-perturbatively. To this aim, meson-meson rescattering is accounted for in a Bethe-Salpeter 
equation which satisfies exact two-body unitarity. 

This method has already been applied successfully to the anomalous decays r]^'^ — > 'y^*^'y^*^ 
and 1]^'^ n'^7i~'y, and to the hadronic decay modes of t] and r]'. The parameters in our 
approach are fixed by the latter two processes and meson-meson scattering phase shifts, so that 
we obtain predictions for the decay widths and spectra of ri,ri' tt'^tt'I'^I' . The decay of r] 
into 7r+7r~e''"e~ is currently under investigation at KLOE@DA$NE and a precise check of our 
prediction for the branching ratio r{rj n~^n~e~^e~)/r{r] — > n~^7r~^) will soon be available [33]. 
Similar investigations are also planned at WASA in Jiilich [2]. 
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A 0{p^) contact term contributions to 77^'^ 7r+7r~7* 

There are several terms in the unnatural parity part of the effective Lagrangian of sixth chiral 
order which contribute to rj^'^ — > 77+77-7* tree level. The full set of Lagrangian terms in the 
SU(3) framework can be found in [34], whereas in the extended U(3) framework — necessary 
to describe 77' decays — the terms relevant for r]^'^ 'K^Ti~'y have been given in [12]. In this 
appendix we repeat the construction of the pertinent Lagrangian terms extending the findings 
of [12] to the description of off-shell photons. 
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The building blocks for the construction of the chiral Lagrangian read 

M ^Wx + X^U , N ^Wx-x'U , (A.l) 

= UW^U , E^, = UWi^DM - (D^DMVU , 

where R^u, are the field strength tensors of the right- and left-handed external fields, 
respectively, the quantity % involves the quark mass matrix Ai — diag(m, m,?^^) (with rh — 
{niu + md)/2), and D^U is the covariant derivative of the meson field U, see [11] for the 
definitions. 

The terms of 0{p^) relevant for the present work are given by 



+ W^{{MC,){C,Q^p) + {N){P,,C^Cp)) 

+ W^{{M{Q^,C^ + C„Qm.)) + {N{P^,C^ - C^P^,))){Cp) 

+ W^o{M){Q,,C^){C0) + Wri{P,,{E\CpC^-C,,C0E\)) 

+ W,,{P,^{E\C,^Cp - CpC^E\)) + W,,{P,,{E\C,^ - C,E\)){Cp) 

+ W^,{P,^{E\Cp-Cf,E\)){C,)} . 

The coefficients Wi are even functions of the singlet field rjo and can be expanded in terms of 
Vo, 

jTrfVo\ _(o) , _(2)»7o , -(4)*7o , /A Q^ 

Wi\^jj ^wl' + wl'— + wyj^ + --- (A.3) 

with expansion coefficients w[''^ not fixed by chiral symmetry. At tree level we find the following 
contribution to the amplitude of rj^'^ — > 7r+7r~7* 

A(^%^'^ - TT+TT-T*) = -ek,e.P^p-,e>^'"''^, /^^g (A.4) 



with 



+ wf^iml - 2ml + 2«+- - k') - w'Si^ml - , 

Pif^ = 327r^y|{8(-4°^ + 3wi'^)ml + (4^°^ + 6^S)(2m^ + ml) 

+ 2^Z;S?(m?' - + 2s+_ - k^) + 3«;S?(mJ + s+_ - k^) 
-2(^S? + 3^iJ)(2m^ -.+_)} . 



(A.5) 
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By defining the combinations 



w^^^ = -2(2«;f ) + + wf^)ml + 8«;f (m^ - ml) + ^ffrnj , 

^i?) = 2wf! + 3«;i? (A.6) 
= -4(2^f - 64°) + wf^ + w?^ + 3^iJ)m' + (4^°^ + 6^S?)(2m^ + ml) , 

4^} = Awf! + 2wf^ + 6^g) + 3^i? , 
= -2^i? - 3^ff , 

which are obviously hnearly independent, we arrive at a simple form for the I3^[,)- 



(A.7) 



Since the mass of the r]' is counted as zeroth chiral order, the w^^^ piece in f3^f^ violates the chiral 

counting scheme. However, as shown in [12], it can be absorbed into the coupling w^^^ and 

in Sec. Owe have employed the renormalized value, f3^f^ = 327r^ 1^2/3(^^7'' +'wl^,'^s+- +w^^'^k'^)^ 
without changing the notation. 
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